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(54) Title: SHARPNESS ENHANCEMENT 

(57) Abstract 

In a method of sharpness en- 
hancement, an input signal (s) is fil- 
tered (HLHP, VLHP) to obtain a fil- 
tered signal (z x , z y ), the filtered signal 
is multiplied (M2, M3) by a control- 
lable fraction (c x , c y ) to obtain multi- 
plied signals, and the multiplied sig- 
nals are added (A2) to the input signal 
(s). According to the invention, the 
controllable fraction (c*, c y ) is gener- 
ated by a non-linear function (HCF, 
VCF). 
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The invention relates to a method and device for sharpness enhancement. 

In scenarios where TV and PC systems merge their functionality, especially in 
the home environment, high-quality displaying of video images in PC architectures becomes a 
challenging issue for industries oriented to the consumer market. Since the spatial resolution 
of synthetic images of the PC outperforms that of natural video scenes, the desire for improved 
quality in TV images will increase. The subjective attribute sharpness, which is determinant by 
the human visual system, is one of the most important factors in the perception of image 
quality. Despite many algorithms for 2-D sharpness enhancement have already been proposed, 
their effectiveness drops in this context. In fact in the image processing literature, it is often 
implicitly assumed that operators that have been designed for the enhancement of 2-D data can 
be straightforwardly applied to video sequences too. As far as contrast sharpening is 
concerned, this is far from true: conventional 2-D techniques introduce many small artifacts 
and non-homogenuities which are perfectly acceptable in a still picture but become very 
visible and annoying in an image sequence. Moreover, the visibility of a defect strongly 
depends on the sequence contents, namely on the amount of details and of motion. 

Moreover, even more attention has to be paid, in the PC environment, to the 
case of decompressed images like MPEG, AVI etc. In such cases, indeed, network and bus 
throughput constraints and limited storage capacity impose the use of compression techniques, 
often operating at relatively low bit rates, which cause visible blocking effects in the 
decompressed images. The artificial high frequencies arising along the borders of the blocks 
themselves are usually amplified by sharpness enhancement algorithms, so that the blocking 
artifact would be emphasized. 

The task of image enhancement techniques is often to emphasize the details of a 
scene so as to make them more visible to a human viewer or to aid some machine performance 
(e.g. object identification); at the same time, they should reduce noise or, at least, avoid its 
amplification. Algorithms for contrast enhancement are often employed in an interactive 
fashion with the choice of the algorithm and the setting of its parameters being dependent on 
the specific application on hand. 
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A large number of approaches have been devised to improve the perceived 
quality of an image [1]. The histogram equalization, a commonly used method, is based on the 
mapping of the input gray levels to achieve a nearly uniform output gray level distribution [1, 
2]. However, histogram equalization applied to the entire image has the disadvantage of the 
5 attenuation of low contrast in the sparsely populated histogram regions. This problem can be 
alleviated by employing local histogram equalization, which is of high computational 
complexity. 

Another method, called statistical differencing, generates the enhanced image 
by dividing each pixel value by its standard deviation estimated inside a specified window 
10 centered at the pixel [2]. Thus, the amplitude of a pixel in the image is increased when it 

differs significantly from its neighbors, while it is decreased otherwise. A generalization of the 
statistical difference methods includes the contributions of the pre-selected first-order and 
second-order moments. 

An alternate approach to contrast enhancement is based on modifying the 
magnitude of the Fourier transform of an image while keeping the phase invariant. The 
transform magnitude is normalized to range between 0 and 1 and raised to a power which is a 
number between zero and one [3]. An inverse transform of the modified spectrum yields the 
enhanced image. This conceptually simple approach in some cases results in unpleasant 
enhanced images with two types of artifacts: enhanced noise and replication of sharp edges. 
20 Moreover, this method is of high computational complexity. 

A simple linear operator that can be used to enhance blurred images is Unsharp 
Masking (UM) [1]. The unsharp masking approach exploits a property of the human visual 
system called the Mach band effect. This property describes the visual phenomenon that the 
difference in the perceived brightness of neighboring region depends on the sharpness of the 
transition [1]; as a result, the image sharpness can be improved by introducing more 
pronounced changes between the image regions. The fundamental idea of UM is to subtract 
from the input signal a low-pass filtered version of the signal itself. The same effect can 
however be obtained by adding to the input signal a processed version of the signal in which 
high-frequency components are enhanced; we shall refer to the latter formulation, 
30 schematically shown in Fig. 1. Fig. 1 shows that a linearly high-pass filtered (HPF) version of 
the input signal x(n) is multiplied by a factor X and thereafter added to the input signal x(n) to 
form the output signal y(n).The output of the high-pass filter HPF introduce an emphasis 
which makes the signal variations more sharp. Its effect is exemplified in Fig. 2. In Fig. 2, 
curve a shows the original signal, curve b shows a first order derivative of the input signal, 
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curve c shows the second order derivative of the input signal, and curve d shows the original 
signal minus k times the second order derivative. 

In one dimension, the UM operation, shown by the block diagram in Fig. 1 , can 
be represented mathematically as: 
5 y(n)=x(n)+Az(n); 

where y(n) and x(n) denote the enhanced signal and the original signal, respectively, z(n) 
means the sharpening component, and X is a positive constant. A commonly used sharpening 
component is the one obtained by a linear high-pass filter which can be, for example, the 
Laplacian operation given by 

10 z(n) = 2x(n) - x(n-l) - x(n+l). 

Even though this method is simple to implement, it suffers from two drawbacks 
that can significantly reduce its benefits. First, the operator introduces an excessive overshoot 
on sharp details. Second, it also enhances the noise and/or digitization effect. The former 
problem comes from the fact that the UM method assigns an emphasis to the high frequency 

15 components of the input, amplifying a part of the spectrum in which the SNR (signal to noise 
ratio) is usually low. On the opposite, wide and abrupt luminance transitions in the input 
image can produce overshoot effects; these are put into further evidence by the human visual 
system through the Mach band effect. 

Several variants of the linear UM technique have been proposed in literature, 

20 trying to reduce the noise amplification. A quite trivial approach consists in substituting a 
bandpass filter for the high-pass one in Fig. 1 . This reduces noise effect, but also precludes 
effective detail enhancement in most images. 

In more sophisticated approaches, Lee and Park [4] suggest to use a modified 
Laplacian in the UM scheme. They propose the order statistic (OS) Laplacian; its output is 

25 proportional to the difference between the local average and the median of the pixel in a 
window. They demonstrate that the resulting filter introduces a much smaller noise 
amplification than the conventional UM filter, with comparable edge-enhancement 
characteristics. A different approach is taken in [5]; in fact they replace the Laplacian filter 
with a very simple operator based on a generalization of the so-called Teager's algorithm. An 

30 example of such an operator is the simple quadratic filter given by 
z(n) =x 2 (n)-x(n-l)x(n+ 1). 

It can be shown that this operator approximates the behavior of a local mean- 
weighted high-pass filter, having reduced high-frequency gain in dark images areas. 
According to Weber's law [1], the sensitivity of the human visual system is higher in dark 
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areas; hence the proposed filter yields a smaller output in darker areas and, therefore, reduce 
the perceivable noise. Even though the above quadratic operators do take into account Weber's 
law, their direct use in unsharp masking may still introduce some visible noise depending on 
the enhancement factor (X in Fig. 1) chosen. 
5 In order to improve the performance of the UM, in [6] the output of the high- 

pass filter is multiplied by a control signal obtained from a quadratic edge sensors: 
z(n) - [x(n-l) - x(n + l)] 2 [2x(n) - x(n-l) - x(n + 1)] (1.1) 
Its purpose is to amplify only local luminance changes due to true image details. The first 
factor on the right-hand side of Eq.1.1 is the edge sensor. It is clear that the output of this 

10 factor will be large only if the difference between x(n-l) and x(n + 1) is large enough, while 
the squaring operation prevents interpreting small luminance variations due to noise as true 
image details. The output of the edge sensor acts as a weight for the signal coming from the 
second factor in Eq.1.1, which is a simple linear high-pass filter. 

Another nonlinear filter, the Rational UM technique, has been devised [7]: the 

15 output of the high-pass filter is multiplied by a rational function of the local input data: 



k[x(n - 1) - x(n + 1)] + h 



In this way, details having low and medium sharpness are enhanced; on the 
other side, noise amplification is very limited and steep edges, which do not need further 
emphasis, remain almost unaffected. Under a computational viewpoint, this operator maintains 
20 almost the same simplicity.8 as the original linear UM method. 

A similar approach is also proposed in [8]; the method is similar to the 
conventional unsharp masking structure, however, the enhancement is allowed only in the 
direction of maximal change and the enhancement parameter is computed as a rational 
function like to the one described in [7]. The operator enhances the true details, limits the 
25 overshoot near sharp edges and attenuates noise in at areas. Moreover it is applied to color 
image enhancement by using an extension of the gradient to multi- valued signal. 

Finally, in [9] the unsharp masking technique is extended with an advanced 
adaptive control that uses the local image content. More precisely, they have concentrated on 
the following properties for adaptive control of the sharpness enhancement: 
30 • local intensity level and related noise visibility; 

• noise level contained by the signal; 

• local sharpness of the input signal; 
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• aliasing prevention, where alias results from nonlinear processing such as clipping. 

These four properties of the video signal are analyzed locally by separate units 
and the amount of the sharpness enhancement depending on this analysis. 

All the proposed contrast enhancement methods have been designed for the 
5 enhancement of 2-D data and cannot be applied to video sequences. In fact they introduce 
many small artifacts and non-homogenuities which are perfectly acceptable in a still picture 
but become very visible and annoying in an image sequence. Moreover, conventional 2-D 
techniques cannot be applied to the enhancement of block-coded image sequences, because 
they emphasize such artifacts that become very visible. The local spectra and bandwidth of 
10 both noise and the signal vary spatially, and the characteristics of the filters need to be locally 
adapted. 



It is, inter alia, an object of the invention to provide an improved sharpness 
enhancement. To this end, the invention provides a shaipness enhancement as defined in the 
independent claims. Advantageous embodiments are defined in the dependent claims. 

In a method of sharpness enhancement according to the invention, an input 
signal is filtered to obtain a filtered signal, the filtered signal is multiplied by a controllable 
fraction to obtain multiplied signals, and the multiplied signals are added to the input signal. 
According to the invention, the controllable fraction is generated by a non-linear function. 

These and other aspects of the invention will be apparent from and elucidated 
with reference to the embodiments described hereinafter. 



In the drawings: 

Fig. 1 shows a linear unsharp masking structure; 

Fig. 2 shows typical behavior of the unsharp masking technique; 

Fig. 3 shows a block diagram of an embodiment of a 3-D unsharp masking 

method; 

Fig. 4 shows a block diagram of an embodiment of a spatial filter for use in the 
embodiment of Fig. 3; 

Fig. 5 shows a block diagram of an embodiment of a control function for use in 
the spatial filter of Fig. 3; 

Fig. 6 shows a plot of a rational function; 

Fig. 7 shows a plot of a temporal controlling function; 

Fig. 8 shows a plot of a control function c x t>; 
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Fig. 9 shows a 3x3 window of a non-linear filter; 

Fig. 10 shows a binary vector that indicates a position of coding artifacts; and 
Fig. 1 1 shows an autocorrelation of the binary vector of Fig. 10. 

5 In this description, we propose a sharpness enhancement technique for video 

applications, taking into account the fact that the images can be blocked. The enhancement is 
accomplished by adding a correction signal to the luminance edges in an Unsharp Masking- 
like way. However, a nonlinear function is used to generate such a correction signal, and both 
spatial and temporal information is used to deal with significant but thin details and with 
10 blocking artifacts; at the same time, a temporal component of the operator reduces the noise 
amplitude. 

This description is organized as follows. In Chapters 2 and 3, a detailed 
description of our new method for video and block-coded image sequences is presented. In 
Chapter 4, we describe an algorithm for the automatic discrimination between video and 

1 5 block-coded image and propose an operator that can be applied for both type of sequences. In 
Chapter 5, our conclusions are given. 

Basically, the enhancement is accomplished by adding overshoot to luminance 
edges in an Unsharp Masking-like way. However, the optimal amount of overshoot added for 
a high image quality depends on the local image statistics. Controls are introduced, to improve 

20 the performances and to adapt it to moving sequences with different characteristics. Both 
spatial and temporal information is also used to enhance details, avoiding blocking artifacts 
and noise. 

2. 3-D unsharp masking technique for video image sequences 

25 

The aim of a contrast enhancement technique is often to emphasize the details 
of a scene so as to make them more visible to a human viewer. Basically, the enhancement is 
accomplished by increasing the steepness of the significant luminance transitions, possibly 
also adding an overshoot to the edges around objects in the image. 
30 In this chapter, we present an Unsharp Masking-based approach for edge 

enhancing for TV application. The proposed scheme enhances the true details, limits the 
overshoot near sharp edges, and attenuates the temporal artifacts. More precisely, to try and 
cope with the problems of linear unsharp masking (noise sensitivity and excessive overshoot 
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on sharp details) we have extended this technique with an adaptive control that uses the local 
spatio-temporal image content. 

The block diagram of the proposed algorithm is depicted in Fig. 3. Let s(n, m, 
t+1) be the input signal; the enhanced signal u(n, m, t) results from the contribution of three 
5 terms: 

u(n, m, t) = s(n, m, t) +Mi s (n, m, t) - u t (n, m, t) (2.1) 
where u s (n, m, t) is the output signal of a spatial non-linear high-pass filter SHPF that 
introduces an emphasis which makes the signal variations more sharp, without incurring 
significant noise amplification and temporal artifacts. u t (n, m, t) is the output of a temporal 

10 high-pass filter THPF operating on three frames s(n, m, t+1), s(n, m, t), s(n, m, t-1), that 

suppresses noise meanwhile preserving the sharpness of moving objects. Delay lines Z" 1 are 
present to generate s(n, m, t) and s(n, m, t-1) from s(n, m, t+1). Multiplier Ml multiplies the 
output signal u s (n, m, t) of a spatial non-linear high-pass filter SHPF by the factor X. Adder A 1 
adds the output signal of the multiplier Ml to the once-delayed input signal s(n, m, t). The 

15 output signal u(n, m, t) is displayed on a display device D. 

The effect of the spatial and the temporal control functions on the image quality 
is explained in Sec.2.1 and Sec.2.1.2 respectively. 

2.1 Spatial operator 

20 The proposed algorithm is formed by a set of directional high-pass filters, 

followed by nonlinear correlation filters. In particular, we consider the separate effects of the 
control signal along the horizontal and the vertical direction: this choice, apart from being 
simple, is also justified by the fact that the eye is more sensitive to the lines and edges having 
these orientations [1], 

25 The block diagram of the spatial high-pass filter SHPF is depicted in Fig. 4, and 

its expression is: 

u s (n, m, t) = z x (n, m, t)c x (n, m, t) +z y (n, m, t)c y (n, m, t) (2:2) 
where 

z x (n, m, t) = 2s(n, m, t) - s(n, m-1, t) - s(n, m + 1, t) (2.3) 
30 z y (n, m, t) = 2s(n, m, t) - s(n-l, m, t) - s(n + l,tn, t) 

respectively, are the outputs of respective Laplacian filters (horizontal and vertical linear high- 
pass filters HLHPF and VLHPF) applied horizontally and vertically to the input image. Their 
amplitude response is a monotone function of the frequency and is used to introduce the 
required emphasis on the signal variation. 
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In order to improve the performance of the high-pass filter we need to condition 
its operation so as to emphasize only the medium-contrast details. To achieve this purpose, z x 
and z y are multiplied by the control functions c x and c y that depend on the local gradient of the 
data. The control functions c x and Cy are obtained by respective horizontal and vertical control 
functions blocks HCF and VCF that each receive the input signal s(n, m, t) of the spatial high- 
pass filter SHPF, z x and z y . Multiplier M2 multiplies the output signal z x (n, m, t) of the 
horizontal linear high-pass filter HLHP by the output c x (n, m, t) of the horizontal control 
function HCF. Multiplier M3 multiplies the output signal z y (n, m, t) of the vertical linear high- 
pass filter VLHP by the output Cy (n, m, t) of the vertical control function VCF. Adder A2 sums 
the outputs of the multipliers M2, M3 to obtain the output signal u s (n, m, t) of the spatial high- 
pass filter SHPF. 



For the sake of simplicity, in the following we explain the operation of our 
method referring only to the x-direction. 



2.1.1 Horizontal control function 

The control function is a nonlinear control that uses the local image content to 
avoid noise amplification, overshoot near sharp edges and temporal artifacts. To better 
understand its behavior, the output of the control function is shown in Fig. 5, which shows the 
20 part of the block diagram enclosed in the dashed box in Fig. 4. 

The control function consists of three major functional blocks: 

• a non-linear rational controlling function RF for the sharpening term; 

• a specialized control TLE for the enhancement of thin lines; 

• an IIR low-pass filtering IIR-LPF of the controlling function. 

25 1x1 the following, the effect of these controls on the image quality will be briefly 

explained. 

Rational function RF 

The controlling function we use is a rational function: 

1 if d x <d x < d. 

30 c x (n,m,t) = \ + dl)d]{n,m ,t) , 



if d x <d t ,d x >d 2 (2.4) 
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The value this function assumes depends on an edge sensor: 
d x (n, m, t) = s(n, m-1, t) - s(n, m+1, t) (2.5) 

It is a band-pass filter that discriminates between signal and noise. To this 
purpose, two implicit hypotheses are made on the significant image details: first, that they are 
5 represented by local gradient values which are larger than those introduced by noise; second, 
that they contribute to the data spectrum mainly in its mid- frequency range. Therefore, this 
function tends to privilege high-gradient areas and it is less sensitive to slow signal variations. 
Moreover, this mechanism makes the proposed operator less sensitive to the Gaussian- 
distributed noise that is always present in the data. 

10 To better understand its behavior, Fig. 6 shows the non-linear control term c x 

as a function of d x for a specific choice of di and d 2 . The main characteristics of the rational 
function c x are: 

1- Cx = 1 for di <d x <d 2 : this characteristic makes the operator able to emphasize 

details that are represented by the low and medium amplitude luminance 
15 transitions; 

2. c x — >Q for d x — ► 0: in homogenous areas, where the noise is usually more 
visible, impulsive noise yields small value of d x (smaller than di), and then a 
small value of the added high-pass signal; 

3. c—Oford* — ► 00: in this way, undesired overshoots near sharp edges are 

20 avoided. In fact, sharp edges will yield high values of d x (larger than d 2 ) and c x 

will be small. 

We can achieve the best balance among these effects by setting the right positions for di and 
d 2 . The intensity of the enhancement can be obviously adjusted by changing the value of X, 

25 Controlling the enhancement of thin lines TLE 

Even though experimental results support the validity of the rational function 
c x , it cannot be used to emphasize thin lines. In fact, to reduce the noise sensitivity, the edge 
sensor is a bandpass filter and is not able to detect lines having a thickness of one pixel. 
Therefore, thin lines are erroneously interpreted as noise and this drawback causes some 

30 artifacts. 

For example, when noisy pixels are adjacent to the line itself, the output of the 
sensor varies along the line. Hence, even for low-amplitude noise, the contribution of the 
correction term c x z x may strongly vary because of the high slope of the control function near 



WO 00/42778 1Q PCT/EPOO/00351 

the origin combined with the very high value of the Laplacian. Temporal changes of the noise 
samples give rise to a non-uniform enhancement, appearing as an annoying flicker effect close 
to the line. 

Thin lines are emphasized in a non-homogeneous way. Moreover, the 
enhancement in the first frame (left) is different from that one in the second frame (right), and 
introduce a temporal artifact which is very visible and annoying in an image sequence. 
Many small artifacts and non-homogenuities are introduced also near lines having slope close 
to, but different from, 0 or 90 degrees. In fact, in digital images, these lines unavoidably show 
jagged borders (staircase effect). In particular, one-pixel-wide lines are formed by segments 
whose width alternates between one and two. The rational function, due to its insensitivity to 
very thin lines, enhances these lines in a non-homogeneous way, and the staircase effect is 
therefore emphasized. 

In order to overcome such problems, a simple but effective method is to include 
a proper control to detect thin lines. To this purpose, we exploit the outputs z x and z y of the 
Laplacian filters themselves. In correspondence of a very thin line, the following three cases 
are possible: 

horizontal line: I z x | = 0 and I Zy | »0 (2 6) 

vertical line: | Z J » 0 and I z y i s 0 (2 .7) 

diagonal line: | Z J » 0 and | z y | » 0 (2 . 8) 

From these relationships, it is straightforward to discriminate between a pixel 
belonging to a line and a noisy pixel, simply by setting a threshold S n which can be derived 
from statistical considerations, as it will be shown in the following. For example, when 

I z x | > S n and I z y | < S„, 
a vertical line is detected and it is only needed to limit the overshoot when such a line is 
already well visible. Therefore, once a thin line is detected by thin line detector TLD, the 
rational function c, is replaced by following overpeaking controlling function OCF (see Fig. 
5) by means of a switch S that is controlled by the output signal of thin line detector TLD: 

177 



Cx = 



which simply limits the correction term z x when its value is larger than d 2 . 

The extension to cases (2.7) and (2.8) is similar. In fact, when the condition 
(2.7) is true, a horizontal line is detected and the vertical rational function c y is replaced by 
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c y . Finally, if condition (2.8) is verified, both c x and c y are replaced by c x and 
c y respectively. 

S n is chosen to discriminate between a pixel belonging to a line and a noisy 
pixel. The value of S n cannot be too small because many noise pixel would be erroneously 
5 interpreted as thin lines; on the other side S n cannot be too much high, otherwise thin line 
would not be individualized. 

A good solution is to determine the relationship between S n and the noise 
variance present in the image. For this purpose, we estimate the probability of erroneously 
labeling a noise pixel as important detail. Supposing a sequence of images degraded with a 
10 zero mean Gaussian noise having variance o 2 , we estimate the probability error when, only 
noise is present. 

Because {s(n + k,m)\k = -l,0,l} and {s{n,m + k)\k = -l,0,l} 

are uncorrected Gaussian random variables with mean zero and variance a 2 (N(0,cr 2 )) , it is 
easy to show that z x and z y are Gaussian random variables, with mean zero and variance 
15 al =a 2 =6<j 2 . 

-x -y 

Since z x and z y depend on the value of s(n, m) (Eq.2.3), they are v.a. correlated. Hence the 
joint probability density function is: 

fzx >z y ( z * >*,) = \ e' Q 



-2/7- 
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g- 1 

2(l-/> 2 ) 
where 

p cov\z x z y \ 

Since Z x ,Z y m N(0,6cr) 

E\Z t Z y \-E\Z x \E\Z y \ 4E\s(n,mf] 4a\ 2 
6a 2 6ct 2 6<t 2 3' 

and the probability of error is: 

P,=\-P<\Z x \{S n ^Z y \{S n ) 

= 1 \ \A^ e -^-"~S' ! >\ d d 
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By selecting S n = 3a ^ we obtain 

P m = 1-0.995 = 0.005= — 
200 

This means that in a 720 x 576 image, for example, more than 2000 noise 
5 points are amplified. However, most of these errors will not be perceived thanks to the 
masking effect of the details and of motion. Experimental results have shown that 

S n = 3a a = 3^6a is a satisfying threshold value. 

We have to observe that this method can be applied only if know the noise 
variance a 2 , but it is not easy to estimate. However, the annoying flicker effects close to the 

10 line and staircase effect are visible above all in those parts of a sequence that are almost 
stationary. In this case, is easy to estimate the noise variance [10]. The simplest method to 
estimate o 2 is frame averaging, which is very effective when the image does not change from 
frame to frame but the degradation does. The simplest and most common form of frame 
averaging is estimating an image f(n, m) from a sequence of N degraded image frame gi(n, m) 

15 fori <i<N. 

Suppose a sequence of degraded images gi(n, m) can be expressed as 
g / (n,m) = /(/i,w) + 7 / («,m), \<i<N (2.9) 
where nj(n, m) is zero-mean stationary white Gaussian noise with variance of a 2 and nj(n, m) is 
independent from ni(n, m) for i^j. If we assume that f(n, m) is nonrandom, the maximum 
20 likelihood (ML) estimate of f(n, m) that maximize 

is given by 

/(»-«) (2-10) 
From (2.9) and (2.10), 

25 /(», m) = /(», m) + -1 £ f/i «) ( 2 - 1 1 ) 

From (2.1 1), the degradation in the frame-averaged image remains a zero-mean 
stationary white Gaussian noise with variance of a 1 I N , which represents a reduction of 
noise variance by a factor of N in comparison with rf i (n 9 m) 
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This suggests estimating 7, from g i . - / . The value of N must be high in order 
to obtain an image f(n, m) with a small noise variance. It cannot be too much high, because 
the estimate of a 2 would be wrong on sequences that are almost stationary. 

Supposing to apply the enhancement operator on sequences with a maximum 
5 value of a 2 =30, we choose N = 5. The estimate of cr 2 will be wrong in conditions of fast 
motion, but in this case the artifacts are less visible. 

IIR low-pass filtering 

From Fig. 5, it can be noticed that the final value of the horizontal control 
10 function c x is obtained by a spatial low-pass filtering, for the reasons that will be explained in 
the following. Formally, its expression is: 

c x (n,mj) = ^c^n-\,mj)+c;(n,mj)] (2.12) 

where c x coincides with c x or c x , depending on whether a thin line is detected or not. 

A non-uniform enhancement effect is present in highly detailed areas with 
15 complex and fine contours, still due to small changes in the output of the sensor amplified by 
the slope of the rational function, together with the high value of the Laplacian in such areas. 
This effect is usually not perceived in 2-D enhancement, thanks to the masking effect of the 
details themselves; however, it generates a flickering which is very annoying in those parts of 
a sequence which are almost stationary. The sharpening is effective, but the edges are 
20 amplified in a non-homogeneous way. 

To make the operator more robust against the appearance of false details, c x is 
IIR filtered as is shown in Eqs.(2.12): the values assumed by on c x neighboring pixels are 
computed by averaging the current value with the previous one along the direction which is 
orthogonal to the control action. In such a way, it is possible to reduce the discontinuity of 
25 these values mainly along the direction of the local borders, avoiding the generation of false 
details and the related visible artifacts. 

It has to be noted that an order of the low-pass filter more than two may 
generate artifacts because in this case we could mediate pixel that don't belong to the same 
object. 

30 



2.1.2 Temporal operator 
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In video sequences, temporal correlation can be exploited to further reduce the 
noise. To suppress noise while preserving the sharpness of moving objects, we propose a 
temporal filter working on three frames. The basic assumption is that when a noisy pixel 
occurs the probability it occurs also in the previous or the next frame in the same spatial 
position is very small. 

Let us define the value of s in three subsequent frames as s(n, m,t-l), 
s(n, m, t) and s(n, m,t + l). When the following conditions 

s(n,m,t-l)== s(n,m,t + l) ^ ^ 

s(n,m,t)- s(n, m, t - 1) * 0 (2 14) 

are both verified, probably a noisy pixel is present in a locally static area. The proposed 
temporal filter is able to evaluate the conditions (2.13) and (2.14) and at the same time to 
perform low-pass filtering. Such a filtering is achieved by subtracting from the input signal the 
following high-pass signal: 

u, («, m, t) = 1 z, (n, m, t)c, («, m, t), (2 j 5) 

1 5 where z, a temporal Laplacian filter: 



10 



20 



25 



30 



z, (w, m,t) = 2s(n, m,t)- s(n, m,t-\)- s{n, m, t + 1), (2 . 1 6) 

and c, controlling function: 

c(^o 4'°"1^ OT ' 0 ^ () '^° 



Its shape is shown in Fig. 7; it can be seen that c, > 0 for \d, \^d 0 ; therefore the temporal low- 
pass filter attenuates signal variations which are smaller than a certain threshold. This allows 
us to smooth selectively the signal variations which are attributed to noise corruption. 

The temporal sensor d, is chosen to discriminate between a noisy pixel and a 
detail. It must be a simple function which is able to evaluate the conditions (2.13) and (2.14). 
To this purpose, the bandpass filter d, =s{n,m,t-\)- s ( n ,m,t + \) gives acceptable results 
with an extremely low computational complexity. As an example, let us consider the case of 
s(n, m, t) corrupted by the noise, no motion, and presence of an uniform area in the masks, so 
that the spatial part of the filter is not allowed to operate: the values of d x , d y and d, will be' 
small, so that c, = c = 0 and c, s 1 
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Hence from Eq.(2.1) the output of the filter becomes: 

* „ t , „ s("> ™, 0 ± s(n, m 9 1 - 1) + s(n 9 m,t + i) 

u{n, m 9 1) = 

4 

which corresponds to a linear low-pass filter that reduces the noise. On the contrary, in 
correspondence of a local movement d t »0 , so that c t = 0 and no filtering is performed. 
5 It has to be noted, however, that this version of the filter may generate artifacts 

in particular conditions of fast motions. Let us consider for example the case of an object 
moving quickly on a uniform background: if the object is present in the filter mask at time t, 
but not at times t-1 and / +7, the value of s(n, m, t) is considered as impulse noise by the 
temporal part of the filter, and the object may be deleted. In an example, a disk with white and 
10 black slices rotates with such a frequency to satisfy conditions (2.13) and (2.14). Applying the 
temporal operator with the bandpass sensor, some artifacts are introduced, which are visible as 
light and dark thin lines. If a larger computational load can be accepted, this type of artifact is 
eliminated by resorting to a Range Filter [11]. Experimental results show that the noise 
reduction is the same with both operators. Therefore, the expression of d t becomes: 
d t {n,m,t) = d m -d mm 
d^ = max{s(/7, m, t - k)\k = -l,0,l} 

d min = min{y(«, m 9 t-k\k = -l,0,l} 

Simulations show that the processed image presents sharper edges which make 
it more pleasant to the human eye. Moreover, the amount of the noise in the image and the 
temporal artifacts are clearly reduced. 

20 

2.2 Computational complexity analysis 

Now, we analyze the computational complexity of our proposed algorithm. 
From Fig. 4 it can be seen that the operator is separable, therefore we consider the separate 
effect of the spatial and temporal operator. The operator needs 19 multiplications and 13 sums. 

25 However, it has to be noted that the values of the non-linear functions c x , c y and c, depend on 
the value of the edge sensors (d x , d y and d t ) which assume only 255 different values. Hence, 
the rational functions can be calculated by using a lookup table; in this case the number of 
multiplications per pixel becomes 7 instead of 19. Hence, under a computational viewpoint, it 
must be stressed that the proposed solution gives acceptable results with an extremely low 

30 computational complexity 
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Unsharp masking tech nique for block-coded image sequences 



3.1 Introduction 

Most video data compression standards use a block-based motion estimate and 
block-based discrete cosine transform (DCT). In particular, they use an N x N pixel-block 
DCT for information packing into a few transform coefficients. This block-based DCT scheme 
takes advantage of the local spatial correlation property of images. However, this type of 
coding induces the well-known blocking artifacts, corner outliers, and ringing noise when the 
image is highly compressed. 

Conventional 2-D techniques for contrast sharpening cannot be applied to the 
enhancement of block-coded image sequences, because they emphasize such artifacts that 
become very visible and annoying in an image sequence. The local spectra and bandwidth of 
both the noise and the signal vary spatially, and the characteristics of the filters need to be 
locally adapted. 

For this purpose, in this chapter we first present a comprehensive 
characterization of the coding artifacts which are introduced into reconstructed video 
sequences through the use of the video coding algorithm and then propose a simple 
modification to the algorithm described in the previous chapter in order to emphasize the 
details and reduces the visibility of the most typical coding artifacts. 

3-2. Video coding distortions 

The following discussion will be limited to descriptions of the artifacts* visual 
manifestations, causes and relationships. We also describe the spatial and temporal 
characteristics of video sequences that are susceptible to each artifact and in which the 
artifacts are visually prominent. 

3.2.1 Blocking effect 

The blocking effect is an annoying side result of techniques in which blocks of 
pixel are treated as single entities and coded separately. In brief, due to the isolated nature in 
which individual blocks are coded, the level and characteristic of the coding error introduced 
into a block may differ from one block to another. So, a smooth change of luminance across a 
border can result in a step in the decoded image. 

When we consider a smoothing operation to remove blocking artifacts in video 
coding, we find three interesting observations. First, the human visual system is more sensitive 
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to blocking artifacts in at regions than in complex regions. Therefore, a strong smoothing filter 
is required on those regions. In complex regions, however, smoothing of few pixel around the 
block boundary is enough to achieve a desired de-blocking effect. Second, smoothing 
operations tend to introduce more undesirable blur in complex regions that in at regions. 
Hence, adaptive smoothing to preserve image details is desirable in complex regions. Third, 
because of motion compensation, blocking artifacts are propagated, and the propagated 
artifacts are more visible in at regions that in complex ones. Therefore, smoothing in that 
regions must cover the inside of a block as well as block boundaries. 

3.2.2 Blurring 

Blurring manifests as a loss of spatial detail and a reduction in sharpness of 
edges in moderate to high spatial activity regions of frames, such as in roughly textured areas 
or around scene object edges. 

For intra-frame coded macroblocks, blurring is directly related to the 
suppressions of the higher-order DCT coefficients through coarse quantization, leaving only 
the lower-order coefficients to represent the contents of a block; therefore, blurring can be 
directly associated with low-pass filtering. 

For predictive coded macroblocks, blurring is mainly a consequence of the use 
of a predicted macroblock with a lack of spatial detail. However, blurring can also be induced 
in bi-directionally predicted macroblocks, where the interpolation of the backward and 
forward predictions results in an averaging of the contents of the final bi-directional 
prediction. In both these cases, the blurred details are supplemented by the prediction error 
that supplies some higher-frequency information to the reconstruction, thereby reducing the 
blurring effect. 

From the observations above, an adaptive enhancement filter is required in 
complex regions in order to emphasize the details and reduce the visibility of the coding 
artifacts. 

3.2.3 Staircase effect 

The DCT basic images are not attuned to the representation of diagonal edges. 
Consequently, more of the higher activity basis images are required to satisfactory represent 
diagonal edges. Due to the typically low magnitude of the higher-order basic images, coarse 
quantization results in their truncation to zero. The contribution originally made by the higher- 
order basic images in forming the diagonal edge is diminished, resulting in a reconstruction 
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exhibiting only the characteristics of the lower frequency basic images, which are generally 
either horizontally on vertically oriented. So, for a block containing a diagonal edge angled 
towards the horizontal, coarse quantization will result in a reconstruction with horizontal 
orientation, and vice versa for blocks angled towards the vertical. Note that the step-wise 
5 discontinuities occur at block boundaries. 



3.2.4 . Ringing 

The ringing effect is related to the high-frequency distortions introduced by the 
DCT coefficient quantization. In fact the representation of a block can be considered as a 

1 0 carefully balanced aggregate of each of the DCT basic images. Therefore, quantization of an 
individual coefficient results in the generation of an error in the contribution made by the 
corresponding basis image to a reconstructed block. Since the higher-frequency basic images 
play a significant role in the representation of an edge, the quantized reconstruction of the 
block will include high-frequency irregularities. This effect is most evident along high contrast 

1 5 edges in areas of generally smooth texture in the reconstruction. 

3.2.5 False edges 
Because of motion compensation, blocking artifacts are propagated; hence, 

false edges are a consequence of the transfer on the block-edge discontinuities formed by the 
blocking effect into the current frame. As with the blocking effect, false edges are mainly 
visible in smooth areas of predictive coded frames. The prediction error in such areas would 
typically be minimal, or quantized to zero, and therefore the false edges would not be masked. 

3.2.6 Motion-compensation mismatch 

25 Motion-compensation (MC) mismatch can be generally defined as the situation 

in which a satisfactory prediction cannot be found for a particular macroblock, resulting in a 
prediction whose spatial characteristics is mismatched with those of the current macroblock. 
The consequence of such a situation is a high level of prediction error which result in a 
reconstruction with highly visible distortions that is a typically a high-frequency noise-like 

30 effect. 
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3-3 Proposed algorithm 

Many de-blocking schemes have been proposed in still image coding such as 
JPEG under the assumption that blocking artifacts are always located at block boundaries. In 
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video coding, however, the blocking artifacts of the previous frame can be propagated to the 
current frame, and can be located at any position in a block due to motion-compensated 
prediction (Sec.3.2.5). Therefore, simple block boundary smoothing is not good enough to 
remove blocking artifacts appearing in video coding. Moreover, a low-pass effect is also 
5 needed inside the macroblock to reduce ringing and motion compensation mismatch effects. 
Iterative methods based on projection on a convex set [12, 13] may be a candidate algorithm. 
However, it is not adequate for real-time video coding due to its complexity. 

In this section, we present an unsharp masking-based approach for noise 
smoothing and edge enhancing. The block diagram is similar to the one shown in Fig. 3 and 
10 Fig. 4, but it is necessary to modify the rational functions (Eq.2.4) in order to enhance true 
details and reduce coding artifacts. 

The new rational function c xb and (c yb ) has to be simple and must obey several 
restrictions, in order to enhances true details and reduce noise: 

1 . It has to attenuate signal variations which are smaller that a certain threshold h. This will 
15 allow us to smooth selectively the signal variations which are attributed to coding artifacts; 

2. The mid-range details which are the signal transitions having magnitude larger than h must 
be enhanced; 

3. To avoid the excessive overshoot over sharp edges: 

20 A simple function verifying all these conditions can be formulated as a rational 

function that is the ratio of two polynomial functions: 



d x (n,m,t) 2 -h 



d x (n,m,t) -h 



7 <1 



c xh {n 9 m 9 t) 



4 kd x (m.,.t) 4 +4Ah' 



kd x (n,m,t) 4 +4 Ah 



otherwise 



d y (n,m,t) 2 -ti 



d y (n 9 m,t) 2 -h* 



4 kd y (m. 9 .t) 4 + 4Ah' 



kd y (n 9 m 9 t) 4 + 4Ah 



otherwise 



where the edge sensors d x and J^are as: 
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a 2 

d x (n,m,t) = d x (n,m,t) 



+<r,h 



„2 (3.3) 

d y (n,m,t) = d y (n,m,t) — - 

CT +< 

Variance-based measurements proved to be a reliable and simple estimator of the degree of 
local activity [14] and have been used for this purpose also in the present work. Hence o 2 
2 is the variance of the pixels belonging to a 3x3 window centered on the current pixel, and 
5 al is a value to be fixed based on the average variance of the image. The parameters it and h ' 
are proper positive factors that are chosen to achieve the best trade-off between enhancement 
of details having low and medium sharpness, reduction of noise and overshoot near already 
sharp edges. The relations between the parameters d,, d 2 . k, h 

1 0 while with the parameter X we can adjust the enhancement intensity. 
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Fig. 8 shows the shape of the horizontal rational control function. It can be seen 
that c xb < 0 for d x <h' ; therefore it attenuates signal variations that are smaller than a certain 
threshold. To better understand the behavior of the filter, we describe the main characteristics 
in at and complex regions. 



1. 
zero ; 



In uniform areas, o 2 is negligible with respect to <r? h , so that d x and d . tend to 
i and c xb ^ c yb < 0 . Hence the output of the spatial filter s + k[z x c xb +z y c yb Jis a low-pass 
signal and the coding artifacts are reduced. 

In particular, supposing that a vertical border occurs between uniform areas 
20 (blocking effect; see Fig. 9); a 2 will be small and c xb - . Since we consider a vertical 
block boundary VB, the previous value of the control function c (n-1 ,m) is similar to . 
Cxb (". m ) ■ Then the output of the IIR filter coincides with c xb (n, m) : 
c x (n,m) = \[c x (n -l,m) + c Xi (k, m)] s c H ( n ,m) = --^. 

Moreover, considering that on such a border z y will be small too, the output of 
25 the filter only depends on the horizontal control function, and we obtain a horizontal linear 
low-pass filter: 
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u s («, m) ~ s(n, m(-\z x (n, m) 

_ 2s(n, m) + s(n — 1, m) + s(w 4- 1, /w) 
= 4 

that reduces the discontinuities between the boundaries of adjacent blocks. 

Our operator also reduces ringing noise. Ringing noise along a vertical edge is 
alike that considered above; in fact the value of z y is small and the output of the filter only 
5 depends on the horizontal control function. So, we obtain a horizontal low-pass filter. 

As regards ringing noise around a diagonal edge, the value of z x and z y are not 
zero and the output of the filter depends on both the horizontal and vertical control functions: 
u s (n, m) =s(n, m) +A.[z x (n, m)c x (n, m) +z y (n, m)c y (n, m)] 

Since a 2 is negligible with respect to <j£ the values of the rational functions are: 

10 C x =C V =— Jr. 

It has to be noted that the previous values of the control functions are greater 
than the current value: 

c x (n -\,m)> c Xh (n, m) and c y (n, m-Y)> c n (n, m) 

because they are closer to the edge. This means that the output of the IIR low-pass filter is 
1 5 greater than - - ! j . Experimental results show that 

c x (n,m) =±[c x (n-\,m) + c Xb (n,m)]= 

c y (n,m) = \ [c y (», m - 1) + c n (n, m)]=-± 

In this way, we obtain a 2-D linear low-pass filter; in fact the expression of u s (n, m) is: 

, . 2s(n, m) + s(n 4- l s m) + s(n - 1, m) + s(n, m + 1) 
u s {n,m)= 

o 

s(n — m — 1) 
+ 6 
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2. In detail areas, o^O, and then d x =d x , d y = d y . Therefore, the mid-range 

details represented by the signal transitions having magnitude larger than h' will be enhanced. 
On the other hand, when d x = 0 {or d y ^ 0) and a 2 » 0 the area is uniform in the horizontal 
(vertical) direction, so that the low-pass filtering is applied. 



25 
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If a smaller computational load is expected, we can calculate the variance of 
five pixels instead of nine pixels. So, if we define 07 like the variance of the following five 
pixels: s(n, m); s(n-l, m); s(n+l, m); s(n, m-1); s(n, m+1), we can use a) in Eq.3.3. 

5 3.4 Conclusions 

From the observations above, the filter for at regions provides a smoothing 
effect inside a block as well as on block boundaries because in video coding the blocking and 
ringing artifacts can be located at any position in a block due to motion-compensated 
prediction. 

0 On the other hand, the behavior of the filter is a nonlinear high-pass one in 

complex areas in order to enhances true details. The variance detector that regulates the filter 
behavior is embedded in the filter itself, which allows a compact expression of the operator. 

It has to be noted that the operator for the enhancement in block-coded image 
sequences is similar to the one described in chapter 2. This suggests defining an enhancement 

5 algorithm for both video and block-coded image sequences. For this purpose, a simple and 
effective solution is proposed in chapter 4. 

i Automatic discriminatio n between video and block-coding image seq uences 

} In chapter 2 we have described an enhancement technique for video image 

sequences that enhances the true details, limits the overshoot near sharp edges and attenuates 
the temporal artifacts. This operator cannot be applied to the enhancement of block-coded 
image sequences, because it emphasizes such coding artifacts that become very visible and 
annoying. 

> For this task, a simple modification has been introduced when processing 
block-coded image sequences in order to emphasize the details and reduce the visibility of 
blocking and ringing effects (see chapter 3). 

However, this operator cannot be applied in video sequences, because it 
introduces significant blurring and does not amplify the low-contrast details. 

> In order to define an enhancement technique for both video and coded image 
sequences, it is important to insert in the operator a mechanism that allows the behavior of the 
filter to change according to the presence of the coding artifacts. To this purpose, in this 
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chapter we insert a weight factor in the expression of edge sensor that depends on a new 
distortion measure for blocking and ringing artifacts. 

4. 1 Enhancement operator for both video and coded image 

From a comparison of the enhancement techniques described in chapter 2 and 3 
it can be seen that they defer for the rational function and the edge sensor. This suggests 
defining a rational function and/or an edge sensor for both video and block-coded image 
sequences. A simple and effective solution is the following: 

d 2 (n,m)-h' d 2 (n,m)-W 



c x (n,m = < 
1 0 where 



kd* (n,m) + 4 Ah 1 kd 4 x («, m) + 4M (4.1) 
1 otherwise 



a 2 

d x («, m) = d x (n, m)— — - (4.2) 

a + P cr ,h 

where the parameter p plays an important role in determining the bias of the filter (the 
criterion for the determination of this weight will be explained in the following). In particular, 
the value of P change according to the presence of the coding artifacts, i.e. P tends to zero in 

15 uncompressed video sequences, while increasing in block-coded sequences. In this way, the 
behavior of the filter becomes more and more low-pass as coding artifacts increases. 

In video sequences, the constant -H in the numerator of c r introduces a low- 
pass effect aiming to reduce noise that is always present in a real image. Moreover, it has to be 
noted that in this case d x -d x , therefore the filter does not introduce significant blurring near 

20 low-contrast details. As mentioned above, the value of P depends on the presence of coding 
artifacts; therefore in the following we introduce a new distortion measure for blocking and 
ringing artifacts. 

4.1.1 Distortion measure for blocking and ringing artifacts 

25 In order to define a metric for blockiness, various methods have been proposed 

in literature [15, 16, 17, 18]. All techniques described suppose that the blocking artifacts are 
always located at block boundaries, but in video coding the blocking artifacts can be located at 
any position in a block due to motion-compensated prediction. Therefore, these methods are 
not good enough to measure the presence of blocking artifacts. A simple solution would be to 

30 apply these algorithms with different value of N; however this approach is not adequate for 
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real-time video coding due to its complexity. Moreover, these techniques do not take into 
account the ringing noise that is very visible and annoying in a image sequences. 
We propose a method to evaluate the presence of both blocking and ringing artifacts that does 
not need information about the position and size of a block. 

5 

Detection of the blocking effect 

. Among the most typical artifacts that can be introduced by coding algorithms, 
the blocking effect is an annoying side result of techniques in which blocks of pixels are 
treated as single entities and coded separately (see Sec.3.3.1). In this case, a smooth change of 
10 luminance across the border of a block can result in a step in the decoded image if neighboring 
samples fall into different quantization intervals. 

Then, if the coefficients of two adjacent blocks are coarsely quantized, we 
expect to see two blocks of low variance and a difference in the intensity slope across the 
block boundary. On the other hand, this abrupt change in intensity slope across the block 
1 5 boundaries of the original unquantized image is rather unlikely. 

From the above consideration, it is clear that an appropriate method to detect a 
blocking effect is to find two adjacent blocks with a very small variance value and a different 
mean value. 

For the sake of simplicity, we consider the horizontal direction. Let s(n, m) be 
20 the luminance level of the pixel in position (n, m); we define , cr^ , /j2 x a\ as the mean 
and variance of two 3x3 adjacent windows: 



(4.4) 



9 



, _ Z ::'.Z£.k+^ J 2 



Z'-'-Z + i,m + 3 + j) 
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2 _ Z Z go [*(" + i,m + +3 + j) - M2t ] 2 



(4.5) 



9 

When the following conditions: 




(4.6) 



are both verified, probably a blocking effect is present in a locally static area. 
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The extension to the vertical direction is trivial. 

Detection of ringing noise 

As mentioned in Sec. 3.2.4, the ringing noise is the Gibb's phenomenon due to 
5 truncation of the high-frequency coefficients by quantization, and, as such, it is most evident 
along high contrast edges in areas generally of smooth texture. It appears as a shimmering or 
rippling outwards from the edge up to the encompassing block's boundary. 

From this consideration, a simple method in order to detect ringing noise is to 
find two adjacent 3x3 windows with a similar mean value (because the ringing noise has 
10 mean zero), but a very different variance value. As done earlier, we consider only the 

horizontal direction; let be , cr* , ju 2r a\ x the mean and variance of two adjacent windows 
(see Eq.4.4). When 



a] =0 

a\ »0 or< 



a\ =0 



&l » 0 (4.7) 



are verified, probably a ringing noise is detected. The ringing noise is most visible as the bit- 
15 rate decrease. 

4.1.2 Determination of p value 

In order to define an enhancement technique for both video and block-coded 
image sequences, we have inserted in the operator a mechanism that allows the behavior of the 
20 filter to change according to the presence of the coding artifacts. 

Equation 4.2 describes the approach that has been chosen in order to achieve 
this goal, where P takes a zero value in uncompressed images, while increasing in block- 
coded images. In this way, the behavior of the filter ranges from a quasi-linear low-pass one in 
uniform area, to a strongly nonlinear high-pass one in detailed areas. Moreover, it becomes 
25 more and more low-pass as coding artifacts increases. 

From the consideration in the previous section it is clear that J3 depends on the 
amount of the coding artifacts in the sequence. To this purpose, we will now show the relation 
between p and the number of coding artifact detected, which take into account the size of the 
image and the number of the frame used for estimate p . 
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For the sake of simplicity, we define t/M vi deo, UM bloC k and UM^ the methods 
described in chapter 2, 3 and the preceding part of chapter 4 respectively. In order to detect 
blocking and ringing artifacts, the conditions (4.6) and (4.7) have to be verified. Experimental 
results show that the following thresholds 

(CT 2 + tT 2 =0 5 
2 ' for detection of blocking effect, and 

\vl-ol. l> 25 

a 2 + a I = | cr 2 - <x 2 2 r | for detection of ringing noise 
are satisfactory. 

In our experiments, we used image sequences after MPEG coding of the 576 x 
720 frames at various bit-rate. Since the blocking effect is more visible in I frames than in B 
10 and P frames due to the motion compensation, the results take into account different frames. 

N b i and N ri indicate how many pixels are interpreted as blocking effect and 
ringing noise respectively. N w = N ri = 0 for uncompressed video sequences, while it increases 
in block-coded sequences as the bit-rate decreases. It has to be noted that, in some sequences 
the ringing noise is more visible than the blocking effect. 
1 5 From the above considerations, it is clear that P depends on the values of N b \ 

and N ri . Obviously, the number of coding artifacts also depends on the size of the image and 
the number of the frame considered; hence let be N c and N r the width and the height of their 
image and N f the number of frame used to calculate N w and N ri , an appropriate measure for p 
is: 

Nbl + Nri (5.1) 

20 B=a K ' 

M NcNrNf 

where a is a constant weight which can be derived from experimental results. In our 
experiment, N c = 720, N r = 576, N f = 4 and a = 3000. 

Our method detects both blocking and ringing artifacts and it does not need any 
25 information about the position and the size of the blocks. Moreover p always take a zero value 
on uncompressed image; therefore is very simple distinguish a coded image between a video 
one. 

In order to further test the noise robustness of the proposed method, an original 
sequence has been corrupted using Gaussian-distributed noise of variance 50. If a prior art 
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linear UM filtering (X = 0.5) is applied, the image is well sharpened, but the noise is amplified 
and it becomes so relevant to degrade of the output scene. The best visual quality is obtained 
using the proposed UM v \dco filtering structure (X = 1.2, di = 30 and d 2 = 40); in fact the spatial 
operator (see Sec.2.1) enhances the medium-contrast details and the temporal filter (see 
5 Sec.2. 1 .2) reduces the noise. 



methods, yields the following results. Using the UM W \^ Q operator, with X = 1, d\= 30 and d 2 = 
40, the processed image is sharp and the noise is almost absent. The C/M b iock method (X = 1, di 
= 30, d 2 = 40, cith 2 = 400) show an improvement of the quality inhomogeneous area; 

10 nevertheless, it introduces significant blurring. In fact, details of intermediate strength are not 
as well defined as in the previous case. Finally, applying the UM^ zpx method to the original 
image, it can be seen that the processed image is similar to that of UM v \dto- hi fact in this 
sequence Nbi = N r i = 0 and f3 - 0. Therefore the low-pass effect of the filter is reduced and the 
low-contrast details are amplified. This means that UM^ m can be applied on both the video 

15 and compressed images. 



Applying a conventional linear UM operator (X = 0.5) to the coded and decoded image, the 

processed image is sharp, but the ringing effect is very visible in the uniform areas. 

The image restored by UM b \ ock operator (X = 1 , di = 30, d 2 = 40, a t h 2 = 400) is free of ringing 

20 effect and the edges are enhanced. The image enhanced by UM b \ ock method is free of blocking 
artifacts. There is practically no blurring of the texture. Further, edges are amplified and the 
image is more sharp. Obviously, the same results are obtained with UM^ apt operator, because 
in compressed images P » 0. However, its effectiveness depends on the value of p. When P = 
0, X = 1, di = 30, di = 40 and a^ 2 = 200, as expected, the sharpening is effective, especially for 

25 medium-contrast details. However, the coding artifacts are amplified in smooth areas. 

Using P = 1, the homogeneous areas are less noisy than in the case of P = 0. Nevertheless, a 
good sharpening is achieved in the detail areas. Finally, with P = 4, it can be seen that coding 
artifacts are removed, but it introduces significant blurring. 



30 has to depend from the presence and the amount of the coding artifacts in the images. Equation 
5.1 describes the approach that has been chosen in order to achieve this goal. 



A comparison of the enhancement between £/Af V ideo. £/M>iock- and £/M a dapt 



The results on compressed image sequences are shown in the following. 



From these considerations, it is clear that P cannot be a constant value, but it 



4.1.3 Observations 
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1 . A major proposed algorithms to define a metric for blockiness suppose that the 
blocking artifacts are always located at block boundaries, but in video coding the blocking 
artifacts can be located at any position due to the motion-compensated prediction. Our 
distortion measure does not need any information about the position and the size of the blocks 

5 and it also detects the ringing noise. 

2. Under a computational viewpoint, it must be stressed that the proposed solution 
maintains almost the same simplicity as the method described in the chapter 3 for the 
enhancement on compressed sequences. In fact, in order to detect blocking and ringing 

10 artifacts we use the variance value which has already been estimated in the expression of the 
rational function (see Eq.3.3). 

3. If our method is applied to the I frames it is possible to determine the size of a 
block; in fact in these frames the blocking and ringing artifacts are always located at block 

15 boundaries. A simple method is described in the following, where for the sake of simplicity, 
we consider only the horizontal direction. First, we define a vector pos x of length N c , where N c 
is the width of the image. When in a block-coded image BCI, a blocking or ringing artifact A 
is detected, the respective position of pos x is set to 1. In this way, a binary vector pos x is 
generated where a 1 value indicates the position of a coding artifact (see Fig. 10). Since the 

20 blocking and ringing artifacts A are located at block boundaries BB, pos x can be expressed as: 



Nc/B 



pos x (n)~£ 5(n-kB) 

/»=! 

where B is the size of the blocks. Then, the autocorrelation of pos x : 
Rx(n)=f; pos x (I) pos x (i+n) 

1=1 

is maximum when n = kB. Therefore, in order to determine the size of the block we can 
25 simply calculate the position of the first maximum of R x and R y . 

As an example, Fig. 1 1 shows the result on a given coded scene after coding at 
1.5 Mb/s. It has to be noted that R x (n) and R y (n) take a zero value for n * k8, therefore the 
image is segmented into 8x8 size blocks. 



30 



4. Our method cannot be applied in graphical image because an abrupt change in 

intensity slope is rather likely; hence many sharp edges would be erroneously interpreted as 



15 
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blocking effect. Moreover, in these images the variance in uniform areas takes a zero value; 
hence the conditions (4.6) is verified for many pixels. 

In order to distinguish a coded image between a graphical one, we observe that 
in the latter the sharp edges can be located at any position; therefore the autocorrelation of 
pos x and pos y takes a high value also for n = 1 (see Eq.4.8): 

Nc 



R x(D « 2 pos x (i) • pos x (i+ 1) 

Nr 

R y0) = X P° s y(0 " pos y (i+l) 



/=0 

From the above consideration, it is simple discriminate the different type of images; in fact the 
following three cases are possible: 
10 (a) N b i + N ri » R x (1) + R y (1) ~ 0: the image is classified as a block coded image; 

(b) N b i + N ri » 0 and R x (l) + R y (1) » 0: the image is classified as a graphical image; 

(c) N b i + N ri ~ 0: the image is classified as a real uncompressed image. 



5. Conclusions 



In this description, we have proposed a sharpness enhancement technique for 
video applications, taking into account the fact that the images can be blocked. 

Conventional 2-D techniques for contrast sharpening cannot be applied to the 
enhancement of video image sequences. In fact, they introduce small artifacts and non- 
20 homogenuities which are perfectly acceptable in a still picture but become very visible and 
annoying in an image sequence; moreover, the visibility of a defect strongly depends on the 
sequence contents, namely on the amount of details and of motion. 

A 3-D UM algorithm has been proposed in chapter 2 in order to overcome such 
problems; despite its simple formulation, it is able to produce a good noise-insensitive 
25 sharpening action without introducing unpleasant overshoot and temporal artifacts. However, 
this operator cannot be applied to the enhancement of block-coded image sequences, because 
it emphasizes such coding artifacts that become very visible and annoying. For this task, a 
simple modification has been introduced when processing block-coded image sequences in 
order to emphasize the details and reduce the visibility of blocking and ringing effects (see 
30 chapter 3). 
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Hence, in coded image sequences, the behavior of the filter ranges from a quasi- 
linear low-pass one in uniform areas, to a strongly nonlinear high-pass one in detailed areas. 
The variance detector that regulates the filter behavior is embedded in the filter itself, which 
allows a compact expression of the operator. 
5 Finally, in order to define an enhancement technique for both the video and 

coded image sequences, a mechanism that allows the behavior of the filter to change according 
to the presence of the coding artifacts, has been inserted in the operator. To this purpose, a 
new distortion measure was introduced which does not need any information about the 
position and the size of the blocks. In fact, the major proposed algorithms to define a metric 
1 0 for blockiness suppose that the blocking artifacts are always located at block boundaries, but 
in video coding the blocking artifacts can be located at any position due to the motion- 
compensated prediction. The algorithm has been evaluated by computer simulations and its 
results clearly show that it outperforms the existing algorithms. A specific component of the 
operator allows this technique to be employed for decoded sequences in a PC environment, 
1 5 also thanks to its low cost in terms of hardware and/or required processing power. The 

algorithm has proven to be robust against the heterogeneous contents of real sequences, giving 
very good performance even on very demanding images. 

It should be noted that the above-mentioned embodiments illustrate rather than 
20 limit the invention, and that those skilled in the art will be able to design many alternative 
embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs placed between parentheses shall not be construed as limiting the claim. The 
word "comprising" does not exclude the presence of elements or steps other than those listed 
in a claim. The word "a" or "an" preceding an element does not exclude the presence of a 
25 plurality of such elements. The invention can be implemented by means of hardware 

comprising several distinct elements, and by means of a suitably programmed computer. In the 
device claim enumerating several means, several of these means can be embodied by one and 
the same item of hardware. The mere fact that certain measures are recited in mutually 
different dependent claims does not indicate that a combination of these measures cannot be 
30 used to advantage. 
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CLAIMS: 



1 . A method of sharpness enhancement, the method comprising the steps of: 

filtering (SHPF) an input signal (s) to obtain a filtered signal (z x , z y ); 
multiplying (Ml, M2, M3) the filtered signal by a controllable fraction (k, c x , 

c y ) to obtain a multiplied signal; and 
5 adding (Al, A2) the multiplied signal to the input signal (s); 

wherein the controllable fraction (c x , c y ) depends on a non-linear function 

(HCF, VCF). 



2. A method as claimed in claim 1, wherein the controllable fraction (c x ) depends 
10 on an edge sensing function (RF). 

3. A method as claimed in claim 2, wherein the controllable fraction (c x ) is smaller 
than its maximum value for edges smaller than a first threshold value (dl). 

15 4. A method as claimed in claim 2, wherein the controllable fraction (c x ) is smaller 

than its maximum value for edges exceeding a second threshold value (d2). 

5. A method as claimed in claim 1 , wherein the non-linear function depends on a 
thin line detection (TLD). 

20 

6. A method as claimed in claim 5, wherein the thin line detection (TLD) is 
carried out on the filtered signal (z x , z y ). 

7. A method as claimed in claim 1, wherein the filtering is a spatial filtering 
25 (SHPF) and the method comprises the further steps of: 

temporally filtering (THPF) the input signal (s) to obtain a temporally filtered 
signal (-u t ); and 

adding (Al) the temporally filtered signal (-u t ) to the input signal (s). 
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8 - A method as claimed in claim 2, wherein the edge sensing function (RF) is 
variance-based. 

9 - A method as claimed in claim 8, wherein the controllable fraction (c x , c y ) is 
5 negative for signal variations than are smaller than a certain threshold. 

10. A method as claimed in claim 1, wherein the method further comprises a 

blocking and/or ringing effect detection including a comparison of mean values of adjacent 
pixel windows. 



10 



1 1 * A method as claimed in claim 1 0, wherein the blocking and/or ringing effect 

detection further includes a comparison of variances of the adjacent pixel windows to a 
threshold value. 



15 1 2 ' A method as claimed in claim 1 0 5 wherein the method further comprises a block 

size detection based on an autocorrelation (R x , R y ) of a vector (pos x , pos y ) indicating blocking 
and/or ringing effects. 



20 



1 3 - A device for sharpness enhancement, the device comprising: 

means for filtering (SHPF) an input signal (s) to obtain a filtered signal (z x , z y ); 
means for multiplying (Ml, M2, M3) the filtered signal by a controllable 
fraction (k, c x , c y ) to obtain a multiplied signal; and 

means for adding (Al, A2) the multiplied signal to the input signal (s); 
wherein the device further comprises means for generating the controllable 
25 fraction (c x , c y ) by means of a non-linear function (HCF, VCF). 



30 



14. 



A television apparatus, comprising: 

a device for sharpness enhancement as claimed in claim 13; and 
means (D) for displaying the sharpened signal. 
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